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ABSTRACT 

In this paper we consider the properties of the inter- 
nal partitions of the nonlinear term, obtained when a filter 
with a sharp cutoff is introduced in wavenumber space. We 
see what appears to be some degree of independence of the 
choice of the position of the cutoff wavenumber for both 
instantaneous and time-integrated partitioned nonlineari- 
ties. We also investigate the basic idea of an eddy-viscosity 
model for subgrid terms and have found that while phase 
modelling will be very poor, amplitude modelling can be 
far more successful. 

INTRODUCTION 

As is well known, full numerical simulation of any sig- 
nificant turbulent flow lies far beyond the scope of cur- 
rent computational resources, the main problem being the 
large number of degrees of freedom involved in the prob- 
lem. As these degrees of freedom may be represented by, 
for instance, the number of independently excited modes in 
wavenumber space, the problem becomes one of eliminat- 
ing modes, in some statistical sense, in order to bring the 
reduced number of degrees of freedom within the capacity 
of current (or even future) computers. One such way by 
which we may systematically obtain such a reduction in the 
number of modes is by the use of a Renormalization Group 
(RG) calculation. A general account of the background to 
this work has been given in the review by McComb (1995). 

In this study, we are undertaking direct numerical sim- 
ulations (DNS) of homogeneous, isotropic, incompressible 
turbulence in a box with periodic boundary conditions, in 
order to assess the underlying feasibility of using RG to 
reduce the size of the computational problem. We have 
already reported some results on the use of conditional av- 
erages (McComb et al. 1997, Machiels 1997) as previously 
formulated by McComb et al (1992) and McComb and 
Watt (1992). In the present paper we concentrate on the 
Hilbert space partitions of the nonlinear terms and their 
filtered projections in order to assess the appropriateness of 
the 'eddy viscosity' concept. Results of this study should 
have direct relevance to large eddy simulations (LES) in 
general, as well as to RG. 



THE PARTITIONED NONLINEAR TERM 

Consider the forced Navier-Stokes equation for station- 
ary turbulence, 

(^■^ + uk^^ Ua(k,t) 

= Mc^^(k) J d^jui3{i,t)u^{'k- i,t) 

+ /c(k,t), (1) 

where u(k, t) is the velocity field in Fourier-space, v is the 
kinematic viscosity, Map-y(k) is given by 

M^l3^{k) = (2i)-' [k^D^^ik) + k.yDc^{k)] , (2) 
where 

D.,(k) = 5.,-^ (3) 

and f(k,t) is a forcing term used to achieve stationarity. 
We may rewrite equation (|l]) in a highly symbolic form as 

Lqu = Muu + f. (4) 

The nonlinear term {Muu in our shorthand notation) 
may be partitioned by introducing a cutoff at fc = fei 
and defining u~ and it'*' such that iia(k, t) = M^(k, t) for 
< A; < fei and Ua{k.,t) = ^^(k, t) for fci < k < ko. The 
maximum cutoff wavenumber, fco, is of the same order of 
magnitude as the Kolmogorov dissipation wavenumber and 
is defined via the dissipation integral, 

e = I uk^E{k)dk ~ / uk^E{k)dk (5) 
Jo Jo 

where e is the dissipation rate. 

Equation (M) can now be expanded to give 

Lou = i,— + f (6) 

where the partitions are defined by 

ip = Mu~u~ (7) 

= 2Mu'u+ (8) 

= Mu+u+ (9) 
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Figure 1: Simulation output. 

(a) skewness ( ), total energy ( ) and dissipation rate 

( ); (b) microscale Reynolds number ( ) and integral 

scale Reynolds number ( ); (c) evolved energy spectrum. 

(te is the eddy turnover time) 
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Figure 2: Correlation between u and at t = 0.5te ( ); 

t = Ite ( ); t = 2te ( ); t = 4ie ( ); t = 6te 

( ); t — 8te ( ) where te is the evolved eddy turnover 

time. The dot-dashed line indicates the position of fc/. 



and we further define 

ip = + + V"^^ = Muu. (10) 

We may now 'solve' equation by introducing Go 
where, for some field X(k,t), 

GoX(k,t) = io'X(k,t) 

-■^''^ii-i')X{k,t')dt' (11) 



so that 



+ 0-+ + 0++ + Go/, 



(12) 



(13) 
(14) 
(15) 



u — 
where 

= GoV-"" 

and we also define 

(j> = (l>~~ + + (?!)++ = Go^A = GoMuu. (16) 

We have carried out direct numerical simulations to cal- 
culate ?/)- and ^-fields in order to investigate their proper- 
ties. 

NUMERICAL SIMULATIONS 

We started from an existing, well validated code for 
the direct numerical simulation (DNS) of turbulence, con- 
structed at the University of Edinburgh and running on the 
Cray T3D administered bv the Edinburgh Parallel Com- 
puting Centre. In Figure hi we have plotted a number of 
the standard DNS outputs generated by our code, running 
on a 256'^ grid. 

For the fundamentals of direct numerical simulations, 
the reader is directed to the pioneering work of Orszag 
(1969 and 1971). Time integration is performed by way 
of a second-order Runge-Kutta scheme and partial dealias- 
ing is achieved through application of a random-shifting 
method (see, for example, Rogallo, 1981). At each time- 
step, the nonlinear term, ip, is calculated by a pseudospec- 
tral method involving a number of fast Fourier transforms. 
ip and may be computed by carrying out the same 
procedure having first zeroed the u"*" and u~ fields respec- 
tively while i/i may be calculated simply by subtracting 
ijj and i/)"*"^ from the total nonlinear term. 



Forcing 

Stationarity is obtained by use of a deterministic forcing 
term given by. 



/.(k,t) 







if < fc < fc/, 
otherwise, 



(17) 



where eo is the desired mean dissipation rate (supplied as 
an input parameter to the simulation), and 



Ef{t) = 



E{k,t)dk 



(18) 



with E{k,t) defined as the energy spectrum, fc/ is chosen 
to be 1.5 so that the forcing is applied to only the first shell 
of wavenumbers. With this forcing, we have observed over 
many simulations that after a sufficient number of time 
steps the velocity field reaches a statistically stationary 
form, as desired. 

Computing the i^-fields 

The evolution of (j) and its partitions, as defined in equa- 
tions ([l3|)-(|l^), is a costly exercise as it must be performed 
in parallel with the evolution of the velocity field. We use a 
simple trapezoidal method to carry out the necessary time 
integrals, but the need to calculate each of the partitions 
of at each time step leads to a code that is roughly three 
times as computationally expensive as a straightforward 
DNS. 

Further problems arise out of the necessity to choose 
suitable initial <ji-fields with which to begin the computa- 
tion. Here, we have chosen to begin by zeroing each of the 
(^-fields for all values of k, in the expectation that after a 
sufficient number of time steps the initial conditions will 
have been forgotten. To this end, we run the DNS code un- 
til the four (ji-fields have reached a statistically stationary 
state, and assume that this indicates convergence to their 
true values. Further evidence may be provided by keeping 
in mind the fact that our forcing is only being applied to 
a single shell in wavenumber space, so that we have 



u — (f) for kf < k < ko- 



(19) 



We note that {|u|'^) « (|</)p) from our evolved data but 
carry out an additional test by computing the correlation 
between the two fields. We define the general correlation 
between two fields, a and b, by 



R{a,b;k) = 



(aa(k)b„(-k)) 
{|a(k)|2)i/2(|b(k)P>V2 



(20) 
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Figure 3: Integrated partitions. 

R{<t>,<^'';k) with fci = 24.5 ( ) and fci = 16.5 (- 

R{(t),(t>' + ;k) with ki = 24.5 ( ) and fci = 16.5 (- 

i?(<^,0++;fc) with fci = 24.5 ( ) and fei = 16.5 (- 

The dot-dashed line indicates k = k\. 



Figure 5: Instantaneous partitions. 

R{ip,ip-';k) with fci = 24.5 ( ) and fci = 16.5 (- 

Rli>,ip- + ;k) with fci = 24.5 ( ) and fci = 16.5 (- 

Rlip,tp++;k) with fci = 24.5 ( ) and fei = 16.5 (- 

The dot-dashed line indicates k — ki. 
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Figure 4: Integrated partitions. 

r{<}),(f)--;k) with fci = 24.5 ( ) and fci = 16.5 (- 

r{4), 4)-+-k) with fei = 24.5 ( ) and fci = 16.5 (- 

r{(j), (j)++;k) with fci = 24.5 ( ) and fci = 16.5 (- 

The dot-dashed line indicates fc = fci. 



Figure 6: Instantaneous partitions. 

r{ip,^--;k) with fci = 24.5 ( ) and fci = 16.5 (- 

rl%lj,^- + -k) with fci = 24.5 ( ) and fci = 16.5 (- 

r(i/),V++;fc) with fci = 24.5 ( ) and fci = 16.5 (- 

The dot-dashed line indicates k — ki. 



and plot R{u, 4>; fc) in Figure ^ for six different time steps. 
We see that by the final time step, the level of correlation 
is excellent for fc > 5 and good for fc > 1. The increasing 
quality of correlation with increasing fc is to be expected as 
a consequence of the fact that higher wavenumbers evolve 
at a greater rate than lower wavenumbers — something 
which is borne out by looking at the correlations computed 
at earlier time steps. The deviation in the first shell is to 
be expected as this is outside the valid range of equation 

Finally we note that in order to compute ^-fields for 
different cutoff wavenumbers, fci, we must reperform the 
entire DNS from initial conditions up to the fully evolved 
state. 

Experimental Details 

We have computed -i/j-fields for a number of different cut- 
off wavenumbers at resolutions of 64^ and 256^ grid points. 
The high cost of calculating the (^-fields means that for 
these we have been restricted to a resolution of 64'^ grid 
points and only two cutoff wavenumbers. 

Our 64'^ simulation achieved a microscale Reynolds num- 



ber of R\ « 70 while our 256^ simulation reached R\ ~ 
190. 

RESULTS 

Low Reynolds Number 

The results given in this section correspond to our 64^ 
simulation. We begin by presenting results for the (f>- 
fields with cutoff wavenumbers fci = 16.5 and fci = 24.5. 
Throughout this work, the cutoff wavenumbers are chosen 
to be half-integers so that they lie between two distinct 
shells. For each data set, we compute the correla 
tween <i> and each of its partitions using equation ( 



ion be- 
\ ) and 

also a measure of their relative magnitudes, r(fc), given by 
(Ib(k)r) 



r(a, b\ fc) 



(21) 



;|a(k)P)- 

Results are plotted in Figures ^ and ^, with both func- 
tions plotted against fc/fci. We first note that, once scaled 
in this way, the exact choice of fci seems to have little effect 
on the shape of the graphs. We believe that the peaky 
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Figure 7: Low Reynolds number (Rx ~ 70) 
R{ip,4' ; fc) for cutoff wavenumbers fci = 4.5, 8.5, 12.5, 
16.5, 20.5, 24.5, 28.5 with ko = 32. The dot-dashed line 
indicates k = ki. 



Figure 10: High Reynolds number {Rx ^ 190) 
/?('(/', for cutoff wavenumbers fci = 16.5, 32.5, 48.5, 

64.5, 80.5, 96.5, 112.5 with ko = 128. The dot-dashed line 
indicates k = ki. 
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Figure 8: Low Reynolds number (Rx ~ 70) 

R{%l),il) '";fe) for cutoff wavenumbers fci = 4.5, 8.5, 12.5, 
16.5, 20.5, 24.5, 28.5 with ko = 32. The dot-dashed line 
indicates k = ki. 



Figure 11: High Reynolds number {Rx ~ 190) 
R{tl',tl'''^;k) for cutoff wavenumbers fci — 16.5, 32.5, 48.5, 
64.5, 80.5, 96.5, 112.5 with ko = 128. The dot-dashed line 
indicates k = ki. 
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Figure 9: Low Reynolds number {Rx ~ 70) 
R{'ip,ip^^ ; k) for cutoff wavenumbers ki = 4.5, 8.5, 12.5, 
16.5, 20.5, 24.5, 28.5 with ko = 32. The dot-dashed line 
indicates k = ki. 



Figure 12: High Reynolds number {Rx ~ 190) 
R{ip,ip^^ k) for cutoff wavenumbers ki = 16.5, 32.5, 48.5, 
64.5, 80.5, 96.5, 112.5 with ko = 128. The dot-dashed line 
indicates k = ki. 



behaviour observed for k/ki < 0.5 in Figure |3| is an effect 
caused by the presence of the forcing term. 

The general picture which seems to emerge from both 
correlation and magnitude information, is that for k <k\, 
the (j>~~ partition is the dominant part of (j). For fe > fei, it 
is tj) ^ which is dominant while (j)^^ is broadly insignificant 
for all values of k. 

It was at this point in our work that it became appar- 
ent that the computational cost involved in calculating the 
(/)-fields was too high and so attention was turned to the 
?/)-fields which are far easier to calculate. For these, a sin- 
gle velocity field realization is enough to calculate ^jj and 
its partitions for any cutoff wavenumber, fci. We begin by 
duplicating the calculations performed on the ^-fields for 
the same resolution grid, and for the same cutoff wavenum- 
bers. Results are shown in Figures ^ and ^ where we see a 
very similar picture to that presented in Figures ^ and |^. 

Taking advantage of the reduction in computational ef- 
fort, we now compute R(;ip, ^ ; k) for a number of differ- 
ent cutoff wavenumbers. The results are shown in Figure [7|. 
We see first of all that there is an excellent collapse of data 
for all cutoffs considered and that there is good correlation 
between ^ and i/) below the cutoff. Above the cutoff, 
this correlation decays rapidly away to zero. We also note 
that, mathematically, (k) = for k > 2fci so that the 
occurrence of a non-zero correlation in this region points 
to the existence of small numerical and aliasing errors. 

This picture is reversed when we consider R{ip,tp~^ ; k) 
as shown in Figure ^ The collapse of data is not so good, 
particularly in the low-wavenumber region, and we see that 
as k/ki increases, the correlation begins to tail away from 
unity. 

A pattern is even more difficult to discern when we com- 
pute R{'ip, ip^'^ ; k) , shown in Figure ^ We see that overall, 
ijj'^^ does not correlate well with tp, although the level of 
correlation increases with k/ki as we move beyond the cut- 
off. 
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Figure 13: , ; fe) for cutoff wavenumbers fci = 

16.5 ( ), 32.5 ( -), 48.5, 64.5, 80.5, 96.5, 112.5 ( ) 

with fco = 128. 

0.005 I ' ^ ' 1 ' 1 

0.004 -I i 
0.003 - ! 




-0.001 - I I 

-0.002 -I I I 

-0.003 - I j - 

-0.004 ]l| i 

-0.005 L ' ' ' ^ ' 

0.00 0.50 1.00 1.50 

k/k, 



Moderate Reynolds Number 

We now extend our work by applying the ideas outlined 
in previous sections to 7/j-field data from our 256^ simula- 
tions. Figures |l(]|-|l2| show correlations of 1/' with each of its 
partitions frornthis data. As can be easily seen, the picture 
has changed very little from our 64"^ data, the biggest dif- 
ference being that the results seem better behaved, which 
is expected due to there being more data points available 
for shell averaging in the region k <k\. 

In the following sections, all results are generated using 
the -i/i-fields taken from our 256'' simulation. 

Eddy-viscosities with Sharp Cutoffs 

We begin by rewriting our Navier-Stokes equation for 
the low-wavenumber modes as 

LoM< =^<--h'</'ir"* + /< (22) 

where the subscript '<' indicates that we are only con- 
cerned with k < ki and where 

1/;^""' = 1/;"+ (23) 

In a large eddy simulation, wavenumbers k > ki will not 

be available and so we introduce some model for 

which we will denote 1/)^ ' . A standard form for ?/)^ ' is 
an eddy-viscosity model, whereby 

ijj^-'^ = -5v{k)k^u< (24) 

for some viscosity increment, 5i'{k). 

We now consider a hypothetical large eddy simulation, 
based around the idea of our wavenumber cutoffs intro- 
duced in previous sections. 



Figure 14: T—{k) ( ) and T^--\k) ( ) for cutoff 

wavenumber k\ — 96.5 with fco — 128. 



We can now compute the correlation between the ex- 
act subgrid terms and the model, R{ip[7~\ip^^~^ ; k). Be- 
cause of -ipl- ^ 's relationship with u — and assuming only 
that 5v{k) is positive for all values of fc — this is equal 

to R{%l>[7~\—u<^;k). The results are plotted in Figure 
|l^ . We see immediately that the correlation is, in general, 
quite poor. We also see what appears to be a difference in 
behaviour between cutoffs fci > 48.5 where there appears 
to be some degree of universality and cutoffs fci < 32.5. 
In general, however, we can state quite categorically that, 
for the case of a sharp cutoff in wavenumber space, no 
eddy-viscosity model (subject to the reasonable constraint 
5v{k) > 0) can perfectly reproduce the missing nonlinear 
terms. But, as we shall see, our conclusion will be modified 
somewhat when we consider phase and amplitude informa- 
tion separately. 

Returning to equation (p^), we can generate an energy 
balance equation by multiplying through by u< and aver- 
aging. Multiplying this through by 27rfc^ then gives, 

^£<-f 2//fc''£;< =r<--f r4""' + w< (25) 

where describes energy transfer to and from the low- 
wavenumber modes through exclusively low-wavenumber 
couplings while T^""' describes energy transfer to and 
from the low-wavenumber modes through coupling involv- 
ing at least one high-wavenumber mode. W<: is the energy 
input due to the forcing. For interest's sake, the two en- 




Figure 15: Su{k) for cutoff wavenumbers fci 
16.5, 32.5, 48.5, 64.5, 80.5, 96.5, 112.5 with ko = 128. 



ergy transfer functions are plotted in Figure We see 
that T has a large negative value in the firstsnell (note 
that the y-axis has been truncated for this graph) corre- 
sponding to transfer of energy away from the energy input 
(forcing). We also see that T is piling up energy at the 
cutoff and that this is balanced by t'^^"' which carries it 
to the higher wavenumbers. 

Introducing an eddy- viscosity model as defined in equa- 
tion (|^ will give us an energy balance equation for LES, 

+ 2pk'^E< + 25v{k)k^E< = T-- + W<. (26) 



dt 



Comparing equations (^5|) and 
form for the eddy viscosity. 



Since we are working in Fourier space, ut will be a com- 
plex scalar and so may be rewritten in the form 

^tT = Ure"'^ (29) 

and similarly for our ■(/'-fields. We may now consider cor- 
relations based solely on phase or amplitude information 
(note that as we are now working with scalar fields, equa- 
tion (|^) must be modified by the removal of the sum over 
components, oK 

In Figures hq — hq we have plotted correlations corre- 
sponding to Figures px|-p^ but concerning ourselves only 
with phase information. NVe see a picture which is largely 
similar to that seen when considering all parts of the %p- 
fields, but note that in places — most obviously for k > ki 
in the -tp'^ data — the correlations are less good. 

Moving on to Figures |l^ — ^ where we have consid- 
ered only amplitude information, we see a different picture. 
Here, while the regions of excellent correlation remain more 
or less untouched, for the rest of the data we see that the 
level of correlation does not drop below about 0.75 and 
indeed seems to remain approximately constant at some 
value between 0.75 and 0.80. 

We can shed some light on the reasons for this by con- 
sidering the statistical nature of the i/i-fields. In Figures 
and 23l we plot probability distribution functions for the 
amphtude and phase components respectively for each of 
the 7/)-fieldqj. We see that the possible values of the ampli- 
tude components are localised and similar in shape (in fact, 
for the chosen value of fc < fci, the PDF's for ipr and tpr~ 
are almost indistinguishable) and hence amplitude correla- 
tions will be good. On the other hand, there is no preferred 
phase and hence phase correlations will be very poor. 

We now carry these ideas across to the analysis of general 
eddy viscosity models outlined in a previous section. Recall 



we can easily derive a the definition of ' as a model for 



/,(— ) 



Siy{k) 



T. 



(--) 



2k'^E^ 



(27) 



For a large eddy simulation, this would ordinarily have to 
be estimated by use of some model, as the whole point of 
LES is the absence of the high-wavenumber modes neces- 



sary for the calculation of T_ 



(--) 



However, with DNS data 



can calculate this, and the results are plotted in Figure 



The general form of these eddy-viscosities appears to be 
in good agreement with the form obtained theoretically by, 
for example, Kraichnan (1976) and we know from previous 
work (see, for example, Lesieur and Rogallo, 1989) and our 
own LES experiments that this particular model provides 
good results. We must now ask why this is, when it is in 
^parent contradiction with the results presented in Figure 

Separating out Phase and Amplitude Effects 

In this section we look at what happens if we separate 
our data into those contributions due to phase and those 
contributions due to amplitude. For each point in our field 
given by the wavevector, k, we choose some unit vector 
perpendicular to k which we call n. We may then generate 
a transverse component of the velocity field. 



Ut 



(28) 



We note that ut provides us with a statistically complete 
picture of our system. This is because, through continuity, 
ul, the component parallel to k will be zero and due to 
isotropy, the statistical properties of ut are independent 
of the exact choice of n. 



■ we can now 

compare the phase and amplitude components of these two 
fields separately. 

In Figure h4 we present the correlation obtained from 
phase information only and see that the general picture 
is the same as when we considered the whole fields. In 
Figure M, however, where we consider only amplitude in- 
formation, we see that there is uniformly good correlation 
at a level of around 0.75. 

This means that the poor correlation between the ex- 
act subgrid terms and an eddy-viscosity model is due al- 
most entirely to mismatching phases. However, by choos- 
ing 5v{k) suitably, it is possible to provide a good match 
as far as amplitude is concerned. 

CONCLUSION 

We have seen that, for the range of cutoff wavenum- 
bers considered here, there is some degree of universal be- 
haviour when we consider the partitions of the nonlinear 
term. We see that for k < ki, the terms involving coupling 
between low-wavenumber modes dominate, and this is a 
useful property when considering the potential for large 
eddy simulations. We have also seen that while it is not 
possible to represent subgrid terms exactly using an eddy- 
viscosity model, it is possible to model the amplitude of 
these terms. This will lead to reasonable results when one 
considers, for example, the fiow of energy but it is clear 
that on a more fundamental level the velocity field will 
become corrupted due to phase-errors. 

Acknowledgements: The results presented in this paper 
have been obtained from simulations performed on the 



^We note that, strictly, these are not PDFs as our data set in 
this instance is insufficiently large — we would expect the true 
PDF of the phase components to be flat, for example. What we 
have instead is a measure of the distribution of that data which 
we do have. 




Figure 16: Phase-correlation. 

R('ipg,ip0~ ; k) for cutoff wavenumbers ki = 16.5, 32.5, 48.5, 
64.5, 80.5, 96.5, 112.5 with ko = 128. The dot-dashed line 
indicates k = k\. 



Figure 19: Amplitude-correlation. 
R{il)r,il'r^ \ k) for cutoff wavenumbers ki = 16.5, 32.5, 48.5, 
64.5, 80.5, 96.5, 112.5 with ko = 128. The dot-dashed line 
indicates k = k\. 
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Figure 17: Phase-correlation. 

R{-tpe,'>pg'^ ; k) for cutoff wavenumbers fci — 16.5, 32.5, 48.5, 
64.5, 80.5, 96.5, 112.5 with ko = 128. The dot-dashed line 
indicates k = k\. 



Figure 20: Amplitude-correlation. 

Ritpr, ipr^'^ k) for cutoff wavenumbers ki = 16.5, 32.5, 48.5, 
64.5, 80.5, 96.5, 112.5 with ko = 128. The dot-dashed line 
indicates k = ki. 




Figure 18: Phase-correlation. 

R{ipe,ipt^;k) for cutoff wavenumbers ki = 16.5, 32.5, 48.5, 
64.5, 80.5, 96.5, 112.5 with ko = 128. The dot-dashed line 
indicates k = ki. 



Figure 21: Amplitude-correlation. 

Riipr, ipr^; k) for cutoff wavenumbers fci = 16.5, 32.5, 48.5, 
64.5, 80.5, 96.5, 112.5 with ko = 128. The dot-dashed line 
indicates k = ki. 
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Figure 22: PDF for Vr- ( ); tpr~ ( ); ( ); 

( ) at A; = 32 with fei = 64.5. 
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Figure 23: PDF for ( ); V-e"" ( ); V-e'^ ( ); 

ip^+ ( ) at A; = 32 with fci = 64.5. 
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Figure 24: Phase-correlation. 

Riii^K ^}o,{i'^<~^}e;k)for cutofF wavenumbers fci = 16.5 

( ), 32.5 ( ), 48.5, 64.5, 80.5, 96.5, 112.5 ( ) with 

ko = 128. 
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Figure 25: Amplitude-correlation. 

ii({V'< ''}r,{V'< ■'}r;) for cutoff wavenumbers ki = 16.5 

( ), 32.5 ( ), 48.5, 64.5, 80.5, 96.5, 112.5 ( ) with 

ko = 128. 



computers of the Edinburgh Parallel Computing Centre. 
The research of A. Young is supported by the Engineering 
and Physical Sciences Research Council. 

REFERENCES 

Kraichnan, R.H., "Eddy viscosity in two and three di- 
mensions," J. Atmos. Sci., Vol. 33, pp. 1521-1536. 

Lesieur, M. and Rogallo, R., 1989, "Large-eddy simu- 
lation of passive scalar diffusion in isotropic turbulence," 
Phys. Fluids A, Vol. 1, No. 4, pp. 718-722. 

McComb, W.D., 1995, "Theory of turbulence", Rep. 
Prog. Phys., Vol. 58, pp. 1117-1206. 

McComb, W. D., Roberts, W. and Watt, A. C, 1992, 
"Conditional-averaging procedure for problems with mode- 
mode coupling", Phys. Rev. A, Vol. 45, pp. 3507-3515. 

McComb, W. D. and Watt, A.G., 1992, "Two-field the- 
ory of incompressible-fluid turbulence", Phys. Rev. A, 
Vol. 46, pp. 4797-4812. 

McComb, W.D., Yang, T.-J., Young, A.J. and Machiels, 
L., 1997, "Investigation of renormalization group methods 
for the numerical simulation of isotropic turbulence," Proc. 
11th Symposium on Turbulent Shear Flows, Grenoble, pp. 
4-23-4-27. 

Machiels, L., 1997, "Predictability of small-scale motion 
in isotropic fluid turbulence," Phys. Rev. Lett., Vol. 79, 
No. 18, pp. 3411-3414. 

Orszag, S., 1969, "Numerical Methods for the Simulation 
of Turbulence," Phys. Fluids (suppl. 2), Vol. 12, pp. 250- 
257. 

Orszag, S., 1971, "Numerical Simulation of Incompress- 
ible Flows Within Simple Boundaries. I. Galerkin (Spec- 
tral) Representations," Stud. Appl. Maths., Vol. 50, No. 
4, pp. 293-327. 

Rogallo, R.S., 1981, "Numerical Experiments in Homo- 
geneous Turbulence," NASA TM-81315. 



